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Abstract 

We study the pairwise annihilation process A+j4 — > inert of a num- 
ber of random walkers, which originally are localized in a small region 
in space. The size of the colony and the typical distance between 
particles increases with time and, consequently, the reaction rate goes 
down. In the long time limit the spatial density profile becomes scale 
invariant. The mean-field approximation of this scenario bears some 
surprises. It predicts an upper critical dimension d c = 2, with log- 
arithmic corrections at the critical dimension and nontrivial scaling 
behavior for d < 2. Based on an exact solution of the one dimensional 
system we conjecture that the mean-field exponents are in fact correct 
even below the upper critical dimension, down to d = 1, while the cor- 
responding scaling function that describes the spatial density profile, 
changes for d < 2 due to the stochastic fluctuations. 

PACS: 82.20.Mj, 05.70.Ln, 02.50.Ey 

Recently there has been much progress made in the understanding of 
the kinetics of diffusing and reacting particles. In particular, the evolution 
of spatial patterns in chemical and biological systems has been subject to 
intensive research for many years. The effects of (chemical) reactions between 
diffusively spreading particles, or of ecological competition on a spreading 
population of species (in the absence of convection) can often be described 
by a reaction diffusion equation of the type 

9 t n = DV 2 n + f(n) (1) 

where n(x, t) is the vector of reactants, D is a matrix of diffusivities, and 
f(n) describes the nonlinear reaction kinetics JlJ. Equation ([]]) is really a 
deterministic approximation ('mean-field'-approximation) of the underlying 
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stochastic process, replacing the discrete and fluctuating occupation numbers 
of the particles by a continuous (coarse grained) population density. In this 
paper we show that for the specific case of a single species of diffusing particles 
A that annihilate each other upon contact A+A —>■ inert, and for a certain set 
of (localized) initial conditions, this mean-field approximation is surprizingly 
powerful: it yields the correct scaling exponents in all dimensions, even down 
to d > 1, while for d < 2 the stochastic fluctuations due to the discreteness 
of the process change only the corresponding scaling function. Generally, in 
contrast, mean field theory is expected to yield exact scaling exponents only 
above the upper critical dimension (which is 2 in this case), wheras below, 
the stochastic fluctuations lead to anomalous exponents. Right at the critical 
dimension we usually find logarithmic corrections to the mean-field behavior. 

For stochastically homogeneous initial conditions, for example, the A + 
A -> inert reaction §, |, §, § 0, g, |, 0, 0, (and the closely related 
reaction A + A — > A, where particles undergo a Brownian motion until they 
hit another particle and undergo coagulation rather than annihiliation, which 
is described by the same mean-field equation as the A + A — > inert process) 
the density decays as n ~ t~ d l 2 for dimensions below the critical dimension 
d c = 2, and as n ~ 1/t for d > 2. The latter formula can be deduced from 
the deterministic ("mean-field") equation of motion for the coarse-grained 
spatial density n(x, t) 

d t n = V 2 n - n 2 , (2) 

where we have absorbed the diffusion constant and the reaction rate into 
the time and length scale. The V 2 -term describes again the diffusion of the 
particles, whereas the n 2 -term accounts for the pairwise annihilation process. 
For a homogeneous density the mean-field equation reduces to n = —n 2 
leading to n ~ 1/t. 

For a certain set of localized initial conditions, however, there is an ex- 
ception to the general change in behavior below the upper critical dimension, 
due to the simple renormalization-group properties of the A + A — > inert- 
reaction pi. In this case the mean- field equation (|2|) itself predicts correct 
exponents, the upper critical dimension and logarithmical corrections: at 
the beginning the system shall be empty up to a finite, populated region. 
From this region particles will spread into the 'vacuum' due to the diffusive 
motion of the particles; and for asymptotically long times we find a scale 
invariant density profile. We will calculate this profile within the mean- field 
approximation and compare it with an exact solution of the one- dimensional 
problem. 

The system under consideration serves as a toy model for systems dis- 
playing a diffusive spreading of reacting agents from a localized region into 
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space, e.g. models for chemical reactions which are ignited at some point. It 
is related to but should not to be confused with the much more complicated 
directed percolation problem [ TTB] , [TJ|, [0| where, starting from a localized 



seed, particles can split, react, die and move diffusively. 

A one-dimensional lattice version of the A + A — > inert-reaction with a 
homogeneous initial state was solved exactly || ||, [10], [□], [12], [16[] owing to 



the fact that the stochastic matrix which enters the master equation can 
be expressed as a bilinear form of fermionic operators. The mean density 
was found to decay as n ~ t _1//2 . In addition, correlation functions were 
derived |9], [TU|, |TT|], confirming the picture of a critical, non equilibrium system 
with anomalous exponents. For dimensions d > 1, an exact RG-approach 
|fj, [7], [Tjj established the upper critical dimension d c = 2 and the power 



law n ~ i" ^ 2 for dimensions d < 2. Consequently, the typical distance 
£ between particles grows like the diffusion length I ~ y 7 ^- This can be 
understood from the recurrence property of random walks for d < 2: Within 
a time intervall t a Brownian particle reacts with particles at distances smaller 



than the diffusion length £, leaving behind voids of linear size ~ £ @, IS , 19 
In order to solve the mean-field equation @ we substitute 

n(x,t) = ^(j)(-^=,tj , (3) 

reflecting the naive dimensions of time, length and the density n. It should 
be kept in mind, that this substitution is without loss of generality. Plugging 
(0) into (0) yields (r = x/\/i is the rescaled position) 

tdttiz, t) = V r 2 + ~ V r • {rcj>) + (1 - d/2) - 2 . (4) 

Introducing a logarithmical time variable s = log(t/to), where to is an arbi- 
trary time scale, we arrive at 

«9 s 0(r, s) = V 2 + ■ (r<j>) + (1 - d/2) - 2 , (5) 

which looks similar to the original equation of motion with, however, a deeply 
modified propagator (linear part). We linearize equation (|) in (f> and split 
it into 

d s <P = n<p+(l-d/2)<p, (6) 

where the operator Q, defined through Q<p = V 2 + (1/2) V • (r0) has zero 
or negative eigenvalues. The unique groundstate (zero mode) of fl is O — 
exp(— r 2 /4); the other eigenfunctions (fi ni ,...,n d = dT ■ ■ ■^d d( Po correspond to 
negative eigenvalues A = — (ni + . . . + n d )/2 (using the notation <9j = d/dri). 
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A critical dimension d c = 2 separates two qualitatively different sectors: 
the linearized equation (||) shows a growing mode for f < d < 2, whereas all 
modes are exponentially damped in s for d > 2. For 1 < d < 2 and a nonzero 
localized initial condition, the unstable mode of will grow until it saturates 
due to the presence of the nonlinear term. There is a unique, stationary 
state 0* of (D, given by V 2 0* + (1/2)V • (r0*) + (1 - d/2)<p* - ((f)*) 2 = 0. A 
closer inspection reveals that 0* is stable, positive, has spherical symmetry 
and decays faster than any power for r — > oo. Going back to the original 
coordinates via @, we obtain the asymptotic time dependence of the density 
n^t) = t~ l 4>* (x/ \ft) . Consequently, in d < 2, the long time decay of the 
total number of particles, for localized initial conditions is given by 

Jd d xn{x,t) ~t d/2 ~\d< 2. (7) 

For d = 2 we see in fact logarithmic terms. The linearized equation (|6]) 
has a marginal mode (zero mode) 0o- Therefore the time dependence of 
its amplitude A is governed by the nonlinear term: d s A oc —A 2 yielding 
A ~ 1/s — l/log(t/t )- Consequently, the total number of particles is not 
constant as a naive extrapolation of (|7|) would suggest, but vanishes according 
to iV ~ l/log(£/£o)- For dimensions above the critical dimension and a 
localized initial condition, behaves like ~ exp((l — d/2)s — r 2 /4) as 
beeing the slowest mode of the linearized equation (^). The nonlinearity 
turns out to be irrelevant since 2 will vanish twice as fast as 0. Using 
the original coordinates, we have n(x,t) ~ t~ d l 2 exp(— x 2 /(4t)), i.e. simple 
diffusion with a conserved total number of particles. In fact, for dimensions 
d < 2 two particles meet with probability one, whereas for dimensions d > 2 
two widely separated particles will miss each other and survive. 

The remainder of this letter is devoted to an exact solution of the stochas- 
tic one-dimensional system along the lines of |5|, |9|, [K| |TT| |12| . We represent 



an occupied site i — . . . , —1, 0, 1, 2, . . . by an occupation number n — 1 and 
an empty site by r« = 0. We start with the initial condition: 

initial state: ... Oil... 1 1 . . . . . 
site: ... —1 1 2 ... L — l L ^ 

where we choose an even number of occupied sites L in order to avoid a 
trivial parity effect. Otherwise, one particle survives - we would see simple 
diffusion of this particle on the long run (This parity effect is absent in the 
mean-field approximation) . 

Now we write down a master equation for the probability P({rj}), as 
done in ||. At first, we define operators i/j and Cj via 

ViP(...,Ti,...) = TiP(. . . ,T h . . .) 
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CjP(. ..,Ti,...) 



= P(...,l- Ti ,...), 



a projection state Q({Ti}) = 1 and an Ising-like operator Sj = 1 — 2z/j yielding 
-1 for an occupied site and +1 for an empty site. The master equation reads 
d t P = M.P with the master operator 

M=Y^ (ciQ+i - 1) (ft + "i+i) • (9) 

i 

It describes the exchange process 01 <-> 10 and the annihilation process 
11 —> 00. The particular ratio between the exchange rate and the annihilation 
rate allows for a quadratic representation in terms of fermionic operators 
but does not constitute a special case since the ratio of both rates is not 
conserved under renormalization j|. We introduce fermionic operators aj = 
. . . •Sj_ 3 Sj_2Sj_iC i i/ i and a\ = . . . •s i _ 3 s i _ 2 s i-i c i(l — ^i) an d find for the master 
operator |5|] 

M = ^2 a \+i a i + <4 a i+i _ 2al a i + 2a i+ iai. (10) 



Following [|T^] we are able to derive the expectation value (z^z/ i+1 ): 

(z^i+i) (*) = (Q\ Wi+i ex P l p o) , (11) 

where P is the initial state (|8|). We have (Q\ Ji4 = and (Q\ q = (Q| and 
therefore 

(fifi+i) (*) = (Q| fli+ifljexp (tA4) |P ) 
= (Q\ exp (-£M) a i+ i exp (tM) exp (-£M) a* exp (tM) \P ) 
= (Q\a l+1 (t) ai (t)\P ), (12) 

where Oj(t) = exp (— t.M) a, exp (t.M). By differentiating this definition with 
respect to t we find easily cij(t) = X),- gi-j(t) a j-, where gi-j(t) = exp(— 2t)Ii-j(2t) 
is the propagator of diffusion on the ID-lattice (Ik is a modified Bessel func- 
tion of order k). Using (|8[) and the definition of cij we arrive at 

L 

(viV i+1 ) = gi+i-j(t)gi- k (t)sign(k - j)(-l) j+k . (13) 
j,k=l 

In order to obtain informations about the particle density we utilize the 
equation of motion for the moments (vj) = ||20|| : 



= (Q\ ViM \Po) = n i+1 + rii-i - 2n< - 2 (viV i+ i) - 2 , (14) 
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where we have used i/jC, = q(1 — z/j) and (Q\ c» = Together with fli~3| ) 
we get 



^ = n i+1 + n 4 _! - 2n; - 2e~ 4t £ (7, +1 _ i (2t) - h-x-^t)) I^ k (2t) 

j,k=i 

x sign(A;- j)(-l) k+j (15) 

and with the help of Jfc_i(2i) — /^ + i(2t) = (k/t)Ik(2t) and after (anti-) sym- 
metrizing j <-> k 

-At L 

dtm = n i+1 + n,_! - 2n< + — - £ ^ - i|/ i _ i (2t)7 i _ fc (2t) (-l) J+fe . (16) 

1 j,k=l 

We are now in a position to perform the continuum limit. We replace the 
discrete position i by a continuous variable x and the discrete Laplacian by 
its continuum version. Furthermore, we approximate the discrete propagator 
of diffusion e~ 2t Ii_j(2t) « (47rt)^ 1/2 exp(-x 2 /(4t)), valid for y/t > L and 
j = 1 . . . L. We have £^=i |Jfe - j|(-l) j+fe = — L for L even and obtain 

d t n{x, t) = d 2 x n(x, t) - exp ■ ( 17 ) 

As before, we find a scale-invariant solution with the ansatz n(x, t) = (p(x/ VI) /t 
(r = x/\/i is again the rescaled position): 

d 2 rV + ^rc^ + if - A exp (-r 2 /2) = 0, (18) 



y,( r ) = L ( e ~ r2 / 2 + - e"^ 4 f dy e"^ 4 ) . (19) 



yielding 

w(r) = 

An V 2 7o 

Figure (0) shows the exact profile p for the one dimensional problem and 
in addition the corresponding mean-field profile, given by an appropriate nu- 
merical solution of the nonlinear differential equation <p"+(r/2)<j)'+<j)—<j) 2 = 0. 
We have adjusted L that both curves have the same integral. Both, the 
stochastic A + A — > inert reaction in ID with an initially localized configura- 
tion of particles and its mean-field approximation do in fact share the same 
scaling exponents. We find scale-invariant profiles n(x, t) = <p*(x/Vt)/t, and 
a total number of particles decaying like N ~ t~ x l 2 . The rescaled profiles 0, 
however, show significant differences. The exact profile has a quartic maxi- 
mum 0(r) oc 1 — r 4 /24, whereas the mean- field profile has a quadratic maxi- 
mum 0(r) oc 1 —0.31 x r 2 . Correlations between different particles, which are 
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Figure 1: Exact rescaled density profile 0(r) (solid), mean-field density profile 
(dashed - both for one spatial dimension). 



neglected in the mean-field approach, are responsible for the difference. For 
dimensions d < 2 strong anticorrelations between particle positions develop 
with time [2^|, since two adjacent particles will recombine with high proba- 



bility. Therefore we have a central region with large voids and a relatively 
low density of particles and a number of particles, which have escaped the 
reaction zone. 

The accuracy of the mean-field approximation generally improves with 
increasing spatial dimension, leading to the conjecture, that the mean-field 
approximation predicts the correct scaling exponents of the spreading colony 
also for dimensions 1 < d < 2 (for dimensions d > 2 mean field is assumed 
to be correct anyway). 

To conclude, we have studied the stochastic A + A — > inert-reaction and 
have demonstrated that the mean-field approximation (0) performs well for 
a particular initial condition. If the particles are located initially in a small 
region in space from which they start to spread out, the mean-field approxi- 
mation apparently predicts the correct scaling exponents of the density profile 
due to the fact, that the renormalization of the A + A — > inert-reaction is 
quite primitive - only the reaction rate but not the propagator is getting 
renormalized j|, |T7| . This should be contrasted to the homogeneous case, 
where mean-field (0) is wrong [23] for d < 2. Of course, in the homoge- 
neous case, the mean-field equation does not know anything about diffusion, 
whereas in the spatially inhomogeneous case the diffusion term is important. 
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Comparing the exact density profile for d = 1 with its mean-field approxima- 
tion, we find significant differences and show that stochastic fluctuations are 
still relevant. It would certainly be interesting to see wether other systems 
with similar RG-properties like the multi-species reaction scheme [f24[] display 
analogous features. 
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